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Abstract 

The present paper is written in response to the comments of M.R. Hadizadeh 
on our original paper ’’Tetraquarks as diquark-antidiquark bound systems” [Phys. 
Lett. B 741 , 124 (2015)]. We present our clarifications on his arguments against 
the accuracy of the procedure and results of the study. The arguments turn out 
to stem mainly from the author’s misunderstanding of the issues discussed in the 
original paper. 


1 Introduction 

In a recent paper titled ’’Tetraquarks as Diquark-Antidiquark Systems” we studied the 
tetraquark systems containing c quark. After revisions and replies to the reviewers’ 
miscellaneous comments, the manuscript was ultimately published in PLB |T]. How¬ 
ever, as research is generative, every discovery brings about further questions. Publi¬ 
cation of ” Tetraquarks as Diquark-Antidiquark Systems” was followed by Comments on 
’’Tetraquarks as Diquark-Antidiquark Systems” by M.R. Hadizadeh [2]. In the hrst place, 
we would like to express our gratitude for his meticulous scrutiny into the article. In the 
second place, we are totally dedicated to presenting our replies to the comments in order 
to compromise on any technical objections as well as misunderstandings apparent within 
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the comments: 


1) As for the coefficient in equations 1, 4, 9 and 11, it must be pointed out that 
it is merely a matter of human mistake while typing. The correct coefficient is ^ . It 
is emphasized that ^ was rightly used in the Fortran code in our calculations, so the 
tetraquark masses were derived based on inserting ^ in the calculations. Indeed, if we 
had used rather than ^ , we would have arrived at different and unacceptable mass 
measures. 


2) The coupling constant we used in Ebert potential has been calculated via (2) in 
Ref. [2]. However, since we made reference to Ebert et ah [3], we didn’t include the 
equation in the body of the article. In our recent work [1], also we employed this coupling 
constant in meson systems. 


3) The following linear variable change is used to transfer the integration interval of 
r/ from [0,rmax] to [-1,+1]: 


r = 


As the author Hadizadeh points out, he has carried out the calculations in conhguration 
space using linear and hyperbolic variable changes, but we conducted our calculations 
using linear variable changes. 

In conhguration space, r-cutoff does not signihcantly inhuence the outcome. The conhn- 
ing term in the potential brings about divergence and thus, we supposed V as constant 
and equal to zero for any r exceeding r-cutoff. In two-body calculations, r-cutoff will have 
no signihcant effects in conhguration space. (In Table 2 of [2], the author Hadizadeh 
presents the fact that varying r-cutoh leads to only minor changes in mass.) However, 
for two- and three-body calculations in momentum space, this renormalization method is 
not applicable, but it works very well with two-body calculations in configuration space 
in long ranges. In our previous works, including [5] (in which M.R. Hadizadeh was an au¬ 
thor), we also successfully employed this method to solve the integral Lippman-Schwinger 
equation and verify the validity of the code and numerical calculations, and we calculated 
the binding energy of ti,bb, and cc systems with high accuracy [5]. Furthermore, besides 
r-cutoh another parameter called r'-cutoh is involved in the calculations of two-body sys¬ 
tems in conhguration space. This parameter can also inhuence the results. 


4) As the author points out, this formalism is applicable to L = 0 states. Some of the 
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details of our methodology in |T] were not given elaborately and this may have triggered 
misunderstanding. As the author holds, we had used a spin-independent unrelativistic 
potential and our results definitely apply to the ground state. Evidently, more accurate 
results are attainable by considering spin effects in the potential. 

However, we found two energies for SA state and three for AA. These were well 
consistent with J 7 ^ 0 state. Hence, we conjecture and suggest that these findings can 
potentially correspond to J 7 ^ 0 state and are in agreement with Eberts [3] results. 

For AA state, there are three angular momentums, i.e. J = and 

J = !■*■“, for states. Thus, utilizing 

Jtot = Jd Jd — {L -\- S)d + (a -|- S’jp 
= 0 ) Jtot = Sd + Sq 

and the parameters M,^X Table 1 of [1], we have calculated different states of J by 
considering the system’s reduced mass. For example, we obtained two energies Ei and 
E 2 {El > E 2 ) yielding A = 1 in for this system’s reduced mass. As a result, Ei is 

supposed to correspond to and E 2 to ■ 

For AA state, we obtained 3 energies Ei,E 2 and E^ {Ei > E 2 > E 3 ) yielding A = 1. 
Thus, we suggest that the highest energy Ei corresponds to J = 2, the lowest energy E 3 
to J = 0, and E 2 to J = 1. 

However, by solving Lippmann-Schwinger equation, we obtained tetraquark masses in 
ground state and propose tetraquark masses for non-zero total angular momentum. As 
stated in the conclusion of pQ, the method enjoys being capable of reducing a complicated 
four-body problem to a simpler two-body problem. Solving the four-body problem via 
the relativistic method for tetraquarks containing light quarks will entail far more com¬ 
plicated calculations. 

Taking all our replies to the disputes against [1] into account, we reject the author’s 
statement that all of the results published in recent letter by M. Monemzadeh et ah [1], 
for masses of tetraquarks with total angular momentum equal to zero, are incorrect and 
the letter has serious numerical problems. Instead, we assert that our results agree with 
other already published results and our method, being simpler, less time-consuming and 
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hence more economical, can be a good replacement for complicated relativistic methods. 
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